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Abstract 

We show that many observable properties of high temperature supercon- 
ductors can be obtained in the frameworks of one-dimensional self-consistent 
model with included superconducting correlations. Analytical solutions for 
spin, charge and superconductivity order parameters are found. The ground 
state of the model at low hole doping is a spin-charge solitonic superstructure. 
Increasing of doping leads to the phase transition to superconducting phase. 
There is a region of doping where superconductivity, spin density wave and 
charged stripe structure coexist. The charge density modulation presents in 
the vicinity of vortices (kinks in the ID model) in the superconducting state. 

PACS numbers: 71.10.Fd, 74.72.-h, 64.60.-i, 71.27.+a 
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Recently discovered stripe phases in doped antifcrromagnets (cuprates and nickelates) [1] 
have attracted attention to the problem of coupled spin and charge order parameters in the 
electron systems. Theoretical [2-5] and experimental [6-10] evidence indicate the possibihty 
that their ground state exhibits spin and charge density waves (SDW and CDW), either com- 
peting, or coexisting with superconductivity. Numerical mean-field calculations [2-4] suggest 
a universality of the spin-charge multi-mode coupling phenomenon in repulsive electronic 
systems of different dimensionahties. Families of the cuprate high-transition-temperature 
superconductors show antiferromagnetism and superconductivity. For the La2-xSrxCu04 
family there is another ordering tendency - unidirectional charge-spin density wave, i.e. 
"stripes". Recent neutron scattering experiment of Lake et al shows that moderate mag- 
netic field makes fluctuating stripes quasi static [6]. An important development in the theory 
of the cuprate superconductors is the prediction that in addition to antiferromagnetism and 
superconductivity there is a tendency toward stripe ordering [2-4] . This prediction is corrob- 
orated by experiments [1,11]- A recent neutron scattering experiment shows that a moderate 
magnetic field can turn a fiuctuating stripe order into a quasi static one in the optimum 
doped cuprates. [7] The vortex state can be regarded as an inhomogeneous mixture of a 
superconducting spin fluid and a material containing a nearly ordered antiferromagnet. 

In this paper we present the one-dimensional effective model describing stripe phase at 
low hole doping and superconductivity state at higher doping. An exact analytical solution 
of the Hartree-Fock problem at and away from half-filling is found. This solution provides 
a unique possibility to investigate analytically the structure of the periodic spin-charge 
solitonic superlattice. It also demonstrates fundamental importance of the higher order 
commensurability effects, which result in special stability points along the axis of concentra- 
tions of the doped holes. Our theory predicts an amazing duality between the spin density 
wave and superconducting order, and implies the presence of stripes near a superconduct- 
ing vortex, and superconductivity near a stripe dislocation. Though there is no long-range 
order in the purely one-dimensional system due to destructive influence of fluctuations, real 
cuprates are three-dimensional, and therefore, the long-range order survives in the ground 
state. Hence, we believe, that one-dimensional mean-field solutions contain universal fea- 
tures of the stripe phase, which are stabilized in higher dimensions. Single-chain analytical 
solutions may be used as building blocks for the stripe and superconducting phase in quasi 
two(three)-dimensional system of parallel chains. 

The Hamiltonian H — Hq + Hg consists of two parts: the Hubbard Hamiltonian with 
on-site repulsion U > 

{i,j)a i i-a 

and the interaction part including superconducting correlations 

i^, = ^A,(04^4 + /..c., (2) 

i 

where is the superconducting order parameter, /j, is the chemical potential. The case of 
the Hubbard model (1) was considered in details earlier [12]. The self-consistent analytical 
solution for the charge-spin solitonic superstructure was found as a function of a hole doping. 
It was shown that effects of commensurability led to a pinning of stripe structure at rational 
filling points \p — 1| = m/n. 
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In the continual self-consistent approximation the effective Hamiltonian can be derived 
similar [12]. We obtain 



A:(-^'-a^+,T + ^+,i^-,T) + 

where A = 2«/7r is a dimcnsionlcss spin coupling constant, As is a dimcnsionless supercon- 
ductor coupling constant, the Pauli matrices, 2a± = dx ± ay, a = U/At; the Plank 
constant is taken as unity, and the length is measured in the units of the lattice (chain) 
period a. In these units momentum and wavevector are dimensionless, and velocity and 
energy posses one and the same dimensionality. The vector ^o-^ = (^o-+, ^o— ) is defined in 
terms of the right- left- moving ^o-±) which constitute the wave function: 

(x) = *+,^e''=^^ + a*_,<,e-''=^^ (4) 

where a = ±1 for a spin | and | respectively. The Fermi-momentum is kp = 7rp/2, 
where in the case of half-filling the average number of electrons per site equals p = 1, i.e. 
kp — 7r/2. The slowly varying real functions A(x) and p{x) are defined as {n{x)) — p{x), 
{S^{x)) = —A{x) cos{nx)/a. The continual approximation requires that a, A, A^ <C 1 (weak 
coupling limit). Note that the constraint A = 2«/7r for the Hubbard model is not necessary 
in a general case, our results remain valid for independent a, A, As <^ 1. 

Introduce p and p as p{x) = p + p{x), J p{x)dx — 0. Then the term ctp^'^^' in Eq. (3) is 
the shift of the chemical potential or the energy, and the term ap^'^'^ can be excluded by 
the unitary transformation (see [12,13]) 



^±{x) — > exp(=piQ; j pdx')^±{x). 



Under this transformation the spin order parameter modifies as A(x) — 

exp(2iQ; pdx')A{x). 

We can diagonalizc the total Hamiltonian H — Hq -\- Hg by performing a unitary Bogol- 
ubov transformations 

'^a{x) = ln,<TUnA^) " <^7n,-a< -al^^) (5) 
n 

which have the form 

n 

in terms of right and left components u±, v± defined as 

Ux) = U,,e''-^ + a/_,,e-^'=-^ f = u,v. (7) 
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New operators 7, 7"*" satisfy the fermionic commutative relations {■jn,a, 'lm,a'} ~ ^m,nSa,a'- 
The transformations (5) must diagonahze the Hamiltonian H: 



en>0 



(8) 



En 



^ en+ [ dx ( 

6„<o 



ttA ttA, 2 , 



where Eg is the ground state energy and e„ > is the energy of excitation n. 
Following [14] we obtain the eigenvalue equations 



Hx = ex, 



(9) 



where 



/ 


-i£ + ap 


A 


A, 





\ 




A* 







A3 






a: 





i|: - ap 


A 




V 





a: 


A* 
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= M_, v+j v.), and self-consistent conditions 

p{x) = 2 + + + v*_v.){l - /)] 



A(x) = -4AEi.l«+/ - - /)] 

A, = 2A,^(l-2/)[t;>+ + 



V u. 



(10) 

(11) 
(12) 



where / = l/(exp[e/T] + 1). We omitted spin indices since in our representation for wave 
functions all equations are diagonal over spin. 

At first, consider homogeneous state with constant A = | A| exp[i(/3], A^ = | A^l cxp[iy?5] 
and p{x) — p = N/L. The average spin density has the form < Sz >oc it!e(A exp(2iA;irx). 
Neglecting trivial dependence on p we obtain two branch spectrum 



with wave functions u,v ca ex^[ikx\ satisfying the symmetry relations 

= ±M_ exp i[ip — (ps]-i V- = ±M+ exp —i[(p + (ps]- 
The self-consistent equations read 

|A| 

where F± = I]J(|A| ± lA^D/e^] tanh[e-i-/T]. At zero temperature we obtain 



(13) 



(14) 



(15) 
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lAI 



log 





-|As| 


|A| 





(16) 



The second equations is derived from the first one by substitution A — * As, A ^ A^. The 
minimum of the ground state energy Eg is achieved at the state A = 2ej?exp[— 1/A], A^ = 
for the case A > A^, and A^ = 2eFexp[— 1/As], A = for the case A < As 

In general case parameters A, As depend on the doping concentration h — \p — 1\. It 
is well known that the coupling constant A monotonically decreases with doping from Aq 
at p = 1 to the value Ao/2 in the limit |p — 1| ^ A/f^ (due to the absence of umklapp 
scattering at p 7^ 1) [15]. If we suppose that superconducting part Hg comes from next 
neighboring site repulsion (as considered for 2D CuO plane model) Hs ~ V pnPn±i-i then 
the self-consistent equation becomes As ~ V n,i^ n±i,}) — 'iV cos kFa{^ i{x)^^{x)) in the 
continual approximation. The coupling constant As ~ |Kcos^ increases with hole doping 
h = 1 — p. If the ground state of undoped system is antiferromagnet state (A > As), phase 
transition to superconducting state will take place at some doping he where A = As. Two 
phases (SDW and SC) with A = As 7^ can coexist at this point. For detailed investigation 
of the phase transition more rigorous consideration of quantum fluctuations is necessary. 

So far we considered uniform state with A, As = const. Since symmetry relations 
between wave function components (14) are independent of absolute values (|A|, |As|), we 
assume that these relations are valid also in a general case of nonuniform order parameters. 
Substituting (14) to (9) we obtain in the case of constant phases (p, (ps equations 



[-icTz-T- + A(7+ + AV_]u 

dx 



eu, 



(17) 



where U"^ = A = (A ± As) exp[i(^]. These equations are eigenvalue equations for 

the Peierls model, were studied in [15,16]. The dependance on p in (17) was excluded 
by means of wave function transformation u±,v± — > exp{^ia J pdx}u±,v±. The term 
a f dxp'^/2 in the total energy Eg is responsible for commensurate effects and pinning of 
the system at rational doping (h = m/n) points [12]. 

Consider the system with A > As. At p = 1 the ground state is antiferromagnet with 
constant A = Aq, p = p and As = 0. As a result of doping kink states are formed 
with local level e = at the center of the gap 2 A. The single kink solution of (17) is 
Ai = A + As = Aotanh(Aox + a/2), A2 = A - A,, = Aotanh(Aox - a/2) with arbitrary 
shift a. The wave functions and the excitation spectrum read 



u± oc (±e -\-k-\- iAotanh^)e 

exp[±i7r/4] 



ifea; ±i7r/4 
c- 1 



= a;2 + a; 



cosh^ ^ 



0, 



(18) 



(19) 



where ^ = AqX ± a/2. The order parameters A, As take form 



An silil 



1 (I 



2(cosh^ Aqx + sinh^ a/2) 



A 



Aq taiili Aoa; 

-j^ _|_ sinh'^ a/2 
cosh"^ A02: 



For the case a = we obtain As = 0, A = Aotanh AqX, p oc 1/ cosh^ AqX. It is a one stripe 
solution found in [12]. The shift < a -C 1 leads to the appearance of the region around 
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the stripe with ^ 0, so that p (x. Xj cosh^ Aqo;, oc a/ cosh^ Aqx. The quasiparticle 
spectrum is independent of a, therefore the equihbrium position a is defined by minimization 
of the potential energy 



TT 



1 1 



Q I ^oQ^te^-\ (20) 
tanh a 4 sinh^ a 

The minimum of the energy (20) is reached at a = for 

For small a the inequality 7 < is possible if A and Ag are not very close to each other 
(I A — As I ^ aAAs). The nontrivial minimum a 7^ exist only in the small region |A — As| ^ 
aWa around the transition point A = A^, where 7 > 0. Stripe and superconductivity phases 
coexist in this region: A^, A,p(a;) 7^ 0. The equilibrium shift a is small a oc ^77 if 7 <^ 1, 
but it 1 logarithmically diverges a oc log 7 in the limit As — > A. 

So we obtain that an increasing of doping for the system with A>As,7<0atp=l 
leads to the forming of the periodic structure of charged kinks A = AotanhAoX, which 
acquires the form at /i = |p — 1| ^ vf/Aq [12] 

A ~ AoV^sn[Aox/v^, k\ p{x) - p oc A^ - A^ 

where K{k) is the Elliptic Integral of the first kind, sn(.,k) is the Jacobi elliptic function, 
\p-l\ = Ao/{2K{k)y/k). 

The superconducting order parameter order Ag 7^ appears in the considered case at 
some higher doping level where 7 becomes positive. In a small region |A — A^j ^ aXXg around 
the transition point A = A^, where 7 > 0, superconductivity and spin/charge orders coexist: 
As,A,p{x) ^ 0. In the particular case of the model with a = this region is reduced to 
the point h — he- A more complicated analysis beyond the scope of the used mean field 
approximation is required at this point to take into account strong quantum fiuctuations, 
including the zero mode due to the degeneration of the ground state with respect to the 
shift a of two sublattices. 

The opposite region A < A^ can be studied using the duality properties of the model. 
It is easy to see that eigenvalues e„ of equations (9) are invariant under transformation 
A < — > Ag. Indeed, if we simultaneously exchange A < — > As and U-{x) < — > v+{x) in Eq. 
(9) the Hamiltonian (9) is not changed (without unimportant terms with p{x)) . Therefore 
the ground state energy Eg in (8) is invariant under the transformations 

A < — > Ag, A < — > Xg. 



Therefore we can apply obtained above solutions for the superconductivity phase. We find 
that in the region 7 > the ground state is superconductor with Ag — const, A = 0. The 
one-dimensional analogue of the vortex in two- or three-dimensional systems is the kink: 
As = AotanhAoo;. 

Due to the duality symmetry the charge density p{x) has the same expression as for the 
kink in spin density wave. Therefore we obtain that the charge density is not zero in the 
vicinity of the kink 
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Similar to the previous case (A > A^j) we find that near the transition point (A ~ A^, 7 > 
0) a stripe structure can arise. In the hmit 0<7-Cl, |p— 1|^ Aq/wf we obtain 

A{x) oc sin7r|p — l\x, p{x) oc cos27r|p — a oc 

To conclude, we have found the self-consistent mean-field analytical solution for the ground 
state structure of the quasi-one-dimensional electronic system with spin, charge and su- 
perconducting correlations. We have found that for an appropriate choice of parameters 
the ground state is striped charge/spin density wave structure at low hole doping. The 
stripe configuration is pinned at rational points |p — 1| = m/n with the pinning energy 
oc exp(— Cn) substantial for small n, which can lead to the stability of the stripe picture. 
The phase transition to the superconductivity state takes place at some doping level. Both 
superconductivity and spin/charge density wave order parameters can coexist in a some 
small region near this point. 

The model is self-dual: The cigcnfunction equations are invariant with respect to trans- 
formations A As, A ^ A^. Therefore properties of superconducting state can be derived 
from the low doping consideration. In particular, we obtained that charge stripes can ex- 
ist as in low doping spin density wave state as in superconducting state in the vicinity of 
spatially nonuniform configurations of A^, for example, vertices (kinks in one dimension). 

Though this one-dimensional model can be applied rather to quasi-one-dimensional sys- 
tems than to high-temperature quasi-two-dimensional anisotropic superconductors, it shows 
some properties peculiar to high-temperature superconductors (one-dimensional stripe struc- 
ture at low doping and superconductivity at a higher doping, etc.). Therefore our results can 
be useful for understanding of high-temperature phenomenon. For describing anisotropic 
properties of real systems a two-dimensional model consideration is required to take into 
account an important contribution from nodal quasiparticles. 

I thank A.V. Balatsky, I. Martin, S.I. Mukhin and J. Zaanen for stimulated discussions 
and K. Machida for a pointing relevant references out. I The work was supported by RFBR 
grant No. 02-02-16354. 
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